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ABSTRACT 

We propose an algorithm for the reconstruction of the signal induced by cosmic strings 
in the cosmic microwave background (CMB), from radio-interferometric data at arcminute 
resolution. Radio interferometry provides incomplete and noisy Fourier measurements of the 
string signal, which exhibits sparse or compressible magnitude of the gradient due to the 
Kaiser-Stebbins (KS) effect. In this context the versatile framework of compressed sensing 
naturally applies for solving the corresponding inverse problem. Our algorithm notably takes 
advantage of a model of the prior statistical distribution of the signal fitted on the basis of 
realistic simulations. Enhanced performance relative to the standard CLEAN algorithm is 
demonstrated by simulated observations under noise conditions including primary and sec- 
ondary CMB anisotropics. 
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1 INTRODUCTION 

Recent observations of the cosmic microwave background (CMB) 
and of the Large Scale Structure of the Universe (LSS) have led 
to the definition of a concordance cosmological model. The anal- 
ysis of the CMB temperature data over the whole celestial sphere 
from the Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite experiment has been playin g a dominant role in designing this 
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According to this model, the cos- 
mic structures and the CMB originate from Gaussian adiabatic 
perturbations seeded in the early phase of inflation of the Uni- 
verse. However, cosmological scenarios motivated by theories for 
the unification of the fundamental interactions predict the exis- 
tence of topologica l defects, resulting from phase transitions at 
the end of inflation JWenkin & Shellardll 19941: lHindmarsh|[T9953 ; 
Iffindmarsh & Kibble! Il995bl : iTurok & Spergeil Il990h . These de- 
fects would have participated to the formation of the cosmic struc- 
tures, also imprinting the CMB. Cosmic strings are a line-like ver- 
sion of such defects which a re also predicted in t he framework 
of fundamental string theory dPavis & Kib"bie1l2005h . As a conse- 
quence, the issue of their existence is a central question in cosmol- 
ogy today. 

Cosmic strings are parametrized by a string tension /i, i.e. a 
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mass per unit length, which sets the overall strength of a string net- 
work. Their main signature in the CMB is characterized by temper- 
ature steps along the string position s. This localized effect, known 
as the Kaiser-Stebbins (KS) effect j Kaiser &Stebbinsl [19841 : iGottl 
hence implies a non-Gaussian imprint of the string network 
in the CMB. The most numerous strings appear at an angular size 
around 1 deg on the celestial sphere. CMB experiments with an 
angular resolution much below 1 deg are thus required in order to 
resolve the width of cosmic strings. 

Experimental constraints have been obtained on a possible 
str ing contribution in term s of upper bounds on the string tension 
u jPerivolaropouloslll993l: iBeyis et alj|2004l ; IWvman et aT]|2005l 
l2006l : lBevis et al.l2007l ; lFraissel2007t) . In this context, even though 
observations largely fit with an origin of the cosmic structures in 
terms of adiabatic perturbations, room is still available for the exis- 
tence of cosmic strings. 

The observed CMB signal can be modeled as a linear su- 
perposition of a statistically isotropic but non-Gaussian string sig- 
nal x proportional to an unknown string tension, with statistically 
isotropic Gaussian noise g comprising the standard components of 
the CMB induced by adiabatic p erturbations, as well as instrumen- 
tal noise. In a first recent work jHammond et alj|2009l) we devel- 
oped an effective method for mapping the string network poten- 
tially imprinted in CMB temperature data, in the perspective of 
forthcoming arcminute resolution experiments performing a map- 
ping of planar patches of the sky in real space, with standard ra- 
diometers. We took a Bayesian approach to this denoising prob- 
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lem, based on statistical models for both the string signal and the 
noise. The denoising is done in the wavelet domain, using a steer- 
able wavelet transform well adapted for representing the strongly 
oriented features present in the string signal. We showed that the 
string signal coefficients are well modeled by generalized Gaus- 
sian distributions (GGD's), which are fitted at each wavelet scale 
using a training simulation borro wed from the set of r ealistic string 
simulations recently produced bv lFraisse et al .1 1 |2008|) . 

In the present work we aim at extending the approach de- 
scribed above to the reconstruction of a string signal from radio- 
interferometric observations of the CMB probing small patches of 
the celestial sphere at the same arcminute resolution. Under the 
most standard assumptions, radio interferometers provide a noisy 
and incomplete Fourier coverage of the planar signal. In these 
terms, the signal reconstruction amounts to solving an ill-posed 
inverse problem encompassing both denoising and deconvolution. 
Relying on the KS effect in the CMB, the string network itself can 
be mapped by the magnitude of the gradient Vi of the string sig- 
nal, which is very sparse or compressible in nature, with significant 
amplitudes only along the strings themselves. 

The theory of compressed sensing offers a new framework for 
solving ill-po sed inverse problem w hen the signa l s are sparse or 
compressible dCandes et al.ll2006d|tJ:ICandedl200"6l ; lDonohdl2006l ; 
iBaraniukl 120071 : iDonoho & Tanned 1 20091) . A band-limited signal 
may be expressed as the N 2 -dimensional vector of its values sam- 
pled at the Nyquist-Shannon rate. By definition, a signal is sparse in 
some basis if its expansion contains only a small number K -C N 2 
of non-zero coefficients. More generally it is compressible if its ex- 
pansion only contains a small number of significant coefficients, 
i.e. if a large number of its coefficients bear a negligible value. The 
theory of compressed sensing demonstrates that, for sparse or com- 
pressible signals, a small number M -C N of random measure- 
ments, in a sensing basis incoherent with the sparsity basis, will 
suffice for an accurate and stable reconstruction of the signals. The 
mutual coherence between two bases may be defined as the max- 
imum complex modulus of the scalar product between unit-norm 
vectors of the two bases. A random selection of Fourier measure- 
ments of a signal sparse in real space or in its gradient are par- 
ticular examples of good sensing procedures. In this framework, 
the reconstruction of the signals may be approached in different 
ways, but all of them encompass a regularization of the ill-posed 
inverse problem by the introduction of a sparsity or compres sibil- 
ity constraint. In a second recent work dWiauxet all l2009ij) . we 
presented results showing that compressed sensing offers power- 
ful image reconstruction techniques for radio-interferometric data, 
which are stable relative to noise, as well as stable relative to non- 
exact sparsity, i.e. compressibility of the signals. These techniques 
are based on global minimization problems, which are solved 
by convex optimization algorithm s dCombettes & Pesquetll2007l ; 
Ivan den Berg & Friedlandenl2008h . We particularly illustrated the 
versatility of the scheme relative to the inclusion of specific prior 
information on the signal in the minimization problems. 

We here wish to apply these results for the reconstruction of a 
str ing signal in the CMB. As for the denoising algorithm proposed 
in lHammond et alj l l2009h . the reconstruction algorithm that we 
present may be considered as a modular component of a larger data 
analysis. It may notably find use as a pre-processing step for other 
meth ods for cosmic string detection based on explicit edge detec- 
tion d-Feong & Smoodl2005l;lLo & Wrightll2005l ; lAmsel et alJI 20071 ; 
iDanos & BrandenbergeJ l2008h . Our reconstruction algorithm is 
tested under different conditions with astrophysical noise compo- 
nents including various contributions to the standard components 



of the CMB, i.e. primary and secondary anisotropics. Our analyses 
rely on the set of re alistic string simulations recently produced by 
iFraisse et ai]d2008T) . 

The remainder of this paper is organized as follows. In Section 
[2] we discuss the string signal and noise in the CMB. In Section[3] 
we describe the inverse problem posed in the context of radio inter- 
ferometry. In Section[4] we define the minimization problems to be 
solved for signal reconstruction in the context of compressed sens- 
ing, notably relying on a prior statistical model of the string signal 
fitted from simulations. In Section [5] we describe our numerical 
analyses and results. We finally conclude in Section|6] 



2 STRING SIGNAL AND NOISE 

In this section, we describe the string signal in the CMB, as well as 
the astrophysical noise made up of the standard primary and sec- 
ondary anisotropics. We also describe the numerical simulations 
used for our subsequent analyses. 



2.1 String signal 

In an inflationary cosmological model, the phase transitions re- 
sponsible for the formation of a cosmic string network occur af- 
ter the end of inflation, so as to produce observable defects. From 
the epoch of last scattering until today, the cosmic string net- 
work continuously imprints the CMB. The so-called scaling solu- 
tion for the string network implies that the most numerous strings 
are imprinted just after last scattering and have a typical angu- 
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strings are also imprinted in the later stages of the Universe evo- 
lution, but in smaller number, according to the number of corre- 
sponding horizon volumes required to fill the sky. 

The main signature of a cosmic string in the CMB is described 
by the KS effect according to which a temperature step is induced 
along the string position. The relative amplitude of this step is given 
by 

ST 
T 

where (3 = w/cand7 = (1— /3 2 ) -1 ' 2 with v standing for the string 
velocity transverse to the line of sight and c for the speed of light, 
and where p is a dimensionless parameter uniquely associated with 
the string tension p through 



(8tt 7 /3) p, 



(1) 



Gp 

n 2 ' 



(2) 



with G standing for the gravitational constant. In the following we 
call p the string tension. 

On small angular scales, realistic simulations can be produced 
by stacking CMB maps induced in different redshift ranges be- 
tween last scattering and today. The simulatio ns that we use in this 
work have been pr oduced by this technique dBouchet et al.|[l988l : 
IFraisse et alj2008h . The string signal is understood as a realization 
of a statistically isotropic but non-Gaussian process on the celestial 
sphere with an overall amplitude rescaled by the string tension p, 
and characterized by a nearly scale-free angular power spectrum: 



C, (p) = p 2 C\ x \ where the positive integer index I stands for 
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Figu re 1. Angular powe r spectra of the string signal and noise (borrowed 
from lFraisse et alj fe008l) ) as functions of the angular frequency in the range 
I € [10 2 ,2 X 10 4 ] in log 10 — log 10 axes scaling. The spectrum of the 
string signal is represented for a string tension p = 2 X 10~ 7 in terms 
of its analytical expression valid at high angular frequencies (red straight 
line). The noise spectra (ordered by decreasing amplitude at low angular 
frequencies) are: the primary CMB anisotropies (black solid line) and the 
gravitational lensing correction (black dashed line), the thermal SZ effect 
in the Rayleigh- Jeans limit (blue solid line), the non-linear kinetic SZ effect 
(green solid line), and the Ostriker-Vishniac effect (magenta solid line). 



the angular frequency index on the sphere. An analytical expres- 
si on of this spec t rum w as provided for I larger than a few hundreds 
bv lFraisse et alj d2008l) . on the basis of their simulations. 

We consider radio-interferometric CMB observations with a 
small field of view corresponding to an angular opening L £ 
[0, 2-7t) on the celestial sphere. In this context, the small portion 
of the celestial sphere accessible is identified to a planar patch of 
size L x L, and we may consider planar signals functions of a 
two-dimensional position vector p. The corresponding spatial fre- 
quencies may be denoted as two-dimensional vectors k with a ra- 
dial component given by the norm k of the vector. In this Euclidean 
limit, the radial component identifies with the angular frequency on 
the celestial sphere, below some band limit B set by the resolution 
of the interferometer under consideration: I = k < B. For k larger 
than a few hundreds, the nearly scale-free planar power spectrum 
of the string signal x(p) reads as: 

P w (fc,p)=p 2 P w (fc), (3) 

with P w (fc) = Cl x) fori = k. 

The observed CMB signal can be understood as a linear su- 
perposition of the string signal and statistically isotropic noise 
of astrophysical and instrumental origin g(p) with some angu- 
lar power spectrum C\ . In the Euclidean limit, the correspond- 
ing planar power spectrum for the noise g(p) may be written as 
p(g) ^ — (jw f Qr i _ f, us n0 [ e th a t we consider zero mean 
signals, identifying perturbations around statistical means. 



2.2 Astrophysical noise 

Current CMB experiments achieve an angular resolution on the ce- 
lestial sphere of the order of 10 arcminutes, corresponding to a limit 
angular frequency not far above B ~ 10 3 . At such resolutions, the 
standard components of the CMB primarily contain the statistically 
isotropic and Gaussian primary anisotropies induced by adiabatic 
perturbations at last scattering. In this context, any possible string 



signal is confined to amplitudes largely dominated by these primary 
anisotropies. The constraints mainly come from a best fit analy- 
sis of the angular power spec trum of the overall CMB signal in 
the WMAP temperature data I Perivolaropo ulos 1993; Bevis et ail 
|2004 IWvman et al.ll2005l, 120061 : iBevis et alj|2007r . lFraissell2007l) . 
The tightest of these constraints l !Fraissel2007l) gives the following 
upper bound at 68 per cent confidence level: 



p p (oxp) = 2.1 x 10" 



(4) 



Forthcoming experiments will provide access to higher angu- 
lar resolution. Among other instruments, radio interferometers such 
as the Arcminute Microkelvin Imager ( AMI) will probe the CMB 
at angular resolutions below 1 arcmin Jjonesll2002l : I Barker et alj 
l2006l : iZwart et al.l2003) . At these resolutions, the secondary CMB 
anisotropies induced by interaction of CMB photons with the 
evolving universe after last scattering will dominate the primary 
anisotropies, and must be accounted for. 

The primary anisotropies exhibit exponential damping at 
high angular frequencies. This contrasts with the slow decay 
of the nearly scale-free angular power spectrum of the string 
signal, which thus dominates over the primary anisotropies at 
high enough angular frequencies. The secondary anisotropies in- 
clude gravitational effects such as the Integrated Sachs-Wolfe 
(ISW) effect, the Rees-Sciama (RS) effect, and gravitational lens- 
ing, as well as re-scattering effects such as the thermal and 
kinetic Sunyaev-Zel'dovich (SZ) e ffects. The SZ effects domi- 
nate these secondary anisotropies dSunvaev &Zel'dovicr]|l98rj 
iKomatsu & Seliakll2002MFraisse et alj2008h . The ISW and RS ef- 
fects associated with the time evolution of the standard gravita- 
tional potentials can be neglected at these angular frequencies. One 
may thus restrict the secondary anisotropies considered in the noise 
to the linear (Ostriker-Vishniac) and non-linear kinetic SZ effects, 
as well as the thermal SZ effect. The SZ effects are actually non- 
Gaussian, spatially dependent, and the kinetic and thermal effects 
are correlated. As a simplifying assumption, we treat these two ef- 
fects as two independent statistically isotropic Gaussian noise com- 
ponents. The effect of gravitational lensing is very small relative to 
the SZ effects, but we still take it into account as a correction to the 
angular power spectrum of the primary anisotropies. 

At arcminute resolution, the thermal and kinetic SZ effects 
have standard deviations around 10 and 5 pK respectively. They 
also have a slow decay at high angular frequencies and will dom- 
inate the string signal for string tension values below the cur- 
rent experimental upper bound. Arcminute CMB experiments are 
in fact primarily dedicated to the detection of these secondary 
anisotropies. Unlike the other effects considered, which have the 
same black body spectrum as the primary anisotropies, the thermal 
SZ effect on the CMB temperature depends on the frequency of 
observation. Its amplitude decreases between the Rayleigh-Jeans 
limit (null frequency) and 217 GHz where it is expected to vanish, 
before increasing again at higher frequencies. Figure [TJ represents 
the angular power spectra as functions of the angular frequency I, 
for a string signal with string tension p — 2 x 10~ 7 , for the primary 
CMB anisotropies and the correction due to gravitational lensing, 
for the Ostriker-Vishniac and the non-linear kinetic SZ effect, and 
for the thermal SZ effect in t he Rayleigh-Jeans li mit. These spectra 
are explicitly borrowed from lFraisse et alj J2008h . 

In this context, the performance of the reconstruction algo- 
rithm to be defined will be studied in the following limits. In a 
first noise condition, we will consider the secondary anisotropies 
as a statistically isotropic Gaussian noise, with power spectrum 
given by the Rayleigh-Jeans limit, that is added to the primary 
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Figure 2. Simulated temperature maps (in /iK) of the string signal in the 
CMB and of its magnitude of gradient for a tension p = 2 X 10~ 7 are rep- 
resented in the top and bottom left panels respectively. Temperature maps 
of the primary CMB anisotropics, i.e. the noise in the PA — IN condition, 
and of its magnitude of gradient are represented in the top and bottom right 
panels respectively. 

anisotropies. In a second noise condition, we will also assume an 
observation frequency around 217 GHz taking advantage of the 
frequency dependence of the thermal SZ effect, and include in 
the noise secondary anisotropies in absence of this effect. This is 
equivalent to including only the kinetic SZ effect and gravitational 
lensing in the secondary anisotropies. Instrumental noise is con- 
sidered to be negligible and simply discarded. These two differ- 
ent noise conditions are respectively denoted as SA + tSZ (sec- 
ondary anisotropies with thermal SZ effect) and SA — tSZ (sec- 
ondary anisotropies without thermal SZ effect) in the following. 
Analyzing these limits can reveal to what extent the kinetic and 
thermal SZ effects hamper the reconstruction of the string signal, 
as a function of the string tension. In a third noise condition, the re- 
construction performance will also be examined in the limit where 
the noise only includes primary anisotropies, neglecting even in- 
strumental noise, and assuming secondary anisotropies have been 
correctly separated. This case is denoted as PA — IN (primary 
anisotropies without instrumental noise) and will be studied in or- 
der to understand the behaviour of the reconstruction algorithm in 
ideal conditions. 

For the sake of our anal yses, foreground emissions such as 
Galactic dust or point sources JKosowskwOOq) are disregarded. 




Figure 3. A map of the primary beam considered is represented in the left 
panel, with darkness increasing with the value of the function A in the range 
[0, 1]. A map of the mask made up of a random selection of 25 per cent of 
the spatial frequencies up to the band limit in the Fourier plane is repre- 
sented in the right panel, with a unit value at the position of the selected 
frequencies identified by black points, and zero otherwise. 



The original maps are sampled at Nyquist-Shannon rate on grids 
with N x N uniformly sampled points pt with 1 i sj A" 2 . We 
consider here samples with N — 128, so that the corresponding 
pixels have an angular size around 0.42', with a band limit B ~ 
2.5 x 10 4 on each each component of the spatial frequencies. 

For each noise component, a simulation may easily be pro- 
duced by taking the Fourier transform of Gaussian white noise, 
renormalizing each Fourier coefficient by the square root of the 
corresponding pow er spectrum, and inverting the Fourier transform 
dRocha et alj|200^) . In each noise condition considered, an over- 
all noise simulation is obtained by simple superposition of the re- 
quired independent components simulated. The power spectrum of 
the noise P' 9 ' (fc) is the sum of the individual spectra. 

Our approach for signal reconstruction notably relies on the 
introduction of specific statistical prior information on the signal in 
the regularization of the ill-posed inverse problem. A number of 64 
simulations of the string signal x are used as training data for fitting 
a prior statistical model for the gradient of the signal. Only 1 addi- 
tional simulation is reserved for testing the algorithm. In all noise 
conditions considered and for each string tension, this test string 
signal simulation is combined with 30 independent realizations of 
the noise in order to produce multiple test simulations. For illustra- 
tion, Figure[2]represents the test simulation of the string signal and 
one simulation of the noise, as well as corresponding maps of the 
magnitude of gradient. 

Note that due to the finite size of the simulations, different 
string signal simulations exhibit slightly different estimated power 
spectra. For simplicity of our argument though both training and 
test simulations have been normalized so as to exhibit the same 
total power for the same string tension (see lHammond et al] ( l2009l) 
for more details). 



2.3 Numerical simulations 

We use simulations of the string signal borrowed fr om the set of 
realistic simulations produced by jFraisse et alj2008h . The specific 
simulations considered are defined on planar patches of size L x L 
for a field of view defined by an angular opening L — 0.9 deg. 
The finite size of the patch induces a discretization of the spatial 
frequencies which may be defined on a grid with N x N uniformly 
sampled points fc; with 1 ^ i ^ A" 2 , below some band limit B. 



3 INTERFEROMETRIC INVERSE PROBLEM 

In this section, we describe the ill-posed inverse problem posed 
for reconstruction of the string signal in the context of radio- 
interferometric observations of the CMB. We also describe a mod- 
ified inverse problem including a whitening operation on the mea- 
sured visibilities and discuss basic ideas for signal reconstruction. 
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3.1 Inverse problem 

Aperture synthesis in radio interferometry is a powerful tech 
radio astronomy dati ng back to more than sixt 



mque in radio astronomy gating bacK to more man sixty 
years ago dRvle & Vonbergill946l: iBlvthd 1 19571: iRvle et alj|l959l : 
iRvle & Hewishlll960l : iThompson et al j|2004> . In this context, the 
small portion of the celestial sphere around the pointing direction 
tracked by a radio telescope array during observation defines the 
original real planar signal or image probed x. The field of view ob- 
served is limited by a primary beam represented by a function A 
with values in the range [0, 1]. Each telescope pair at one instant 
of observation identifies a baseline defined as the relative position 
between the two telescopes. To each baseline is associated one mea- 
surement called visibility. In the simplest setting one also considers 
baselines with negligible component in the pointing direction of the 
instrument. Under this additional assumption, if the signal is made 
up of incoherent sources, each visibility corresponds to the value 
of the Fourier transform of the signal multiplied by the primary 
beam at a spatial frequency k identified by the components of the 
baseline in the plane of the signal. 

We assume that the spatial frequencies k probed by all tele- 
scope pairs d uring the observation belong to the discrete grid of 
points ki (see lWiaux et alj J2009ah for more details). The Fourier 
coverage provided by the M/2 spatial frequencies probed kb, with 
1^6^ M/2, can simply be identified by a binary mask in 
the Fourier plane equal to 1 for each spatial frequency probed and 
otherwise. The visibilities measured may be denoted by a vec- 
tor of M/2 complex Fourier coefficients y £ C M//2 = {yb = 
y(kb)}i^b^M/2- Formally, the visibilities may equivalently be de- 
noted by a vector of M real measures y £ R M = {y r }i^r^M 
consisting of the real and imaginary parts of the complex measures. 
In this discrete setting, the Fourier coverage is in general incom- 
plete in the sense that the number of real constraints M is smaller 
than the number of unknowns JV 2 : M < iV 2 . An ill-posed inverse 
problem is thus defined for the reconstruction of the signal, seen as 
a vector of real sampled values x £ M. N = {x(pi)} 1 ^ i ^ N 2 , from 
the measured visibilities y. In the context of our analysis, the noise 
affecting the signal consists in the vector of real sampled values for 
the standard CMB components g £ ~R N = {g(pi)}i^i^N 2 ■ The 
inverse problem therefore reads as: 



<t>x + n with <t> = MFA 
and with n = <t>g, 



(5) 



where the measurement matrix <t> £ q(m/2)xN l( \ en iif[ es the 
complete linear relation between the signal and the visibilities. The 
matrix A £ K Jv2xJv2 = {Aij = A(pi)5ij} liitj<N 2 is the di- 
agonal matrix implementing the primary beam. The unitary matrix 



F £ 



i/N} 



N 2 implements the 



discrete Fourier transform providing the Fourier coefficients. The 
matrix M £ R( M / 2 ) xiv2 = {MyJ^^Mp.^ji is the rect- 
angular binary matrix implementing the mask characterizing the 
interferometer. It contains only one non-zero value on each line, 
at the index of the Fourier coefficient corresponding to each of the 
spatial frequencies probed kt- Neglecting small correlations intro- 
duced between spatial frequencies by the primary beam, the power 
spectrum of the noise g is nonetheless modified as P^ A9 \k) = 



\A(k)\ 2 P<- 3 \k). The 



C M ' 2 = {n,,}, 



s^b^M/2 



af- 



fecting the visibilities is Gaussian noise with different variances 
p( Ag \kb) for different components rib- Let us recall that g is 
real so that the components of n are independent only when as- 
suming that the mask M contains no pair of antipodal spatial fre- 



quencies {kb, —kb}. For each such pair, the two corresponding 
noise components are complex conjugate of one another. In this 
work M/2 denotes the number of independent complex visibili- 
ties, while Mo/2, with Mo ^ M, will stand for the total number 
of complex visibilities including possible pairs of antipodal spatial 
frequencies. 



3.2 Whitening 

The measurement matrix may be augmented by a whitening op- 
eration which will be useful in a reconstruction perspective. It 
consists in dividing each measured visibility yb by the standard 
deviation of the real and imaginary parts of the corresponding 
noise component rib- This operation is simply implemented in 
terms of the diagonal matrix W £ r(w/2)x(m/2) = = 

[P (Ag) (kb)/2] 1/2 8 b b'}i^b,b'^(M/2)- This gives rise to the modi- 
fied inverse problem 

y = i>x + n with <J> = WMFA 

and with h = <t>g, (6) 

where the whitened visibilities y £ C M//2 = {yb = 
y(kb)}i^b^M/2 are now affected by identically distributed inde- 
pendent (i.i.d.) Gaussian noise with unit variance of the real and 
imaginary part of each component rib, i- e - white Gaussian noise. 

For each test simulation, the visibilities are computed in ac- 
cordance with relation QJ. The primary beam A considered for our 
observational set up is a Gaussian function with a full width at half 
maximum (FWHM) equal to the angular opening of the field of 
view FWHM — L. The function has a maximum value equal to 
unity. Two Fourier coverages are considered independently, defined 
by masks made up of a ra ndom selection of Mo spatial frequencies 
(see lWiauxetall j2009ij) for more details) with either 25 per cent 
or 50 per cent of the Fourier plane probed up to the band limit 
B. Note that, after discarding redundant visibilities associated with 
pairs of antipodal spatial frequencies in the two masks considered, 
these two percentages respectively correspond to selections of M 
spatial frequencies with 21.7 and 37.4 per cent. The primary beam 
as well as the mask for a Fourier coverage of 25 per cent are illus- 
trated in Figure [3] The visibilities are also whitened in accordance 
with relation l[6j by the introduction of the suitable matrix W in 
each noise condition. 



3.3 Reconstruction basics 

In the perspective of the reconstruction of the signal x, relation 
l(6j represents the measurement constraint. Considering a candi- 
date reconstruction x, the residual noise reads as n £ 

c M/2 = 

{^i,}i^i,^A//2 — V — ^x. The residual noise level estimator, de- 
fined as twice the negative logarithm of the likelihood associated 
with x, simply reads as 

«/2 L-(*s)J 2 



X 



(y,<t>,x) = 2Y / 



P( A s) (k b ) 



(7) 



The notation | ■ | for a scalar stands for the complex modulus when 
applied to a complex number and for the absolute value when ap- 
plied to a real number. The £2 norm of the residual noise is the stan- 
dard norm of the corresponding vector: 1 1 Tl 1 1 2 = (Y2b=l l^b] 2 ) 1 ^ 2 - 
This noise level estimator follows a chi-square distribution with M 
degrees of freedom. Typically, this estimator should be minimized 
by the good candidate reconstruction. The measurement constraint 
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on the reconstruction may be defined as a bound x 2 (y, t t > , a^) ^ e 2 , 
with e 2 corresponding to some suitable large percentile of the x 2 
distribution in order to avoid noise over-fitting. 

In this context, many signals may formally satisfy the mea- 
surement constraint. A regularization scheme that encompasses 
enough prior information on the original signal is needed in order 
to find a unique solution. All image reconstruction algorithms will 
differ through the kind of regularization considered. 

Firstly, the simplest possible reconstruction procedure would 
consist in assuming that the non-probed spatial frequencies have a 
null value, and simply inverting the Fourier transform after dividing 
each visibility by the variance of the noise in order to enhance the 
signal-to-noise ratio. This reconstruction can by written as xq oc 
( WM F) t y, where the notation stands for the conjugate transpose 
operation. This map simply identifies with the convolution of the 
original signal and noise A(x + g) by a beam (F t M t W 2 MF)(5o, 
where So stands for a spike at the center of the field of view: do £ 
R^ = {^ojji^j^Ara, with io identifying the central pixel. In 
the standard vocabulary of radio interferometry, the beam is called 
the dirty beam, and the convolved map is called the dirty map in 
recognition of its visual appearance contrasting with the original 
signal. A normalization constant needs to be chosen so as to get 
a unit central value of the beam. In other words, no denoising or 
deconvolution is performed by this simple procedure that we call 
INVERT. The only important operation in this reconstruction is 
the enhancement of the signal-to-noise ratio by the matrix W 2 . 

Secondly, the standard CLE AN algorithm is essentially a 
Matc hing Pursuit (MP) algorithm dMallat & Zhandll993l ; TMallatl 
Il998h which performs the deconvolution by iterative removal of 
the dirty beam, thus leading from the dirty map to a clean map 
jHogbom|[l974l : ISchwad[l978l : iThompson et alj|2004l) . Formally, 
the algorithm should be stopped when the \ 2 constraint above is 
satisfied. 




Figure 4. Positive tail of modeled prior GGD distribution V(\7x\p) for a 
string tension p = 2 X 10 — 7 . The GGD is superimposed on the histogram 
of the aggregate of the magnitude of the gradient of the 64 test simulations. 



the real signal under a \ constraint on the candidate reconstruction 
x: 



mm 



|x||tv subject to x (y, <t>, x) ^ e 2 



(8) 



The vector of sample values of the magnitude of the gradient 
associated with the vector x is denoted as Va; 6 R^ = 
{Vx(pi)} 1 ^ is g JV 2 . By definition, the TV norm of a signal is the 
£i norm of the magnitude of its gradient, i.e. ||:e||tv = ||Va;||i = 

I Va~(pi) I- We set t0 be m e 99 th percentile of the \ 2 dis- 
tribution. 

The iterative algorithm considered for solving this convex 
optimiza tion problem is based on the Douglas-Rachford s plitting 
method dCombettes & Pesquedl2007l : iFadili & Starckl|2009h in the 
framework of proximal operator theory jMoreadl 1962 ). 



4 COMPRESSED SENSING RECONSTRUCTION 

In this section, we describe the standard global minimization prob- 
lem TVe for the reconstruction of the string signal in the context of 
compressed sensing. We model the prior statistical distribution of 
the gradient of the string signal with a GGD, fitted on the basis of 
the training simulations. We also propose to enhance the TV £ prob- 
lem on the basis of this statistical prior information, hence defining 
a STV e minimization problem. We finally propose a power spec- 
tral model (PSM) used for independent estimation of the string ten- 
sion, which may notably be useful for accelerating the convergence 
of the iterative algorithms solving the minimization problems. 



4.1 TV e problem 

For the reconstruction of a signal of sparsity K, the theory of com- 
pressed sensing proves that the number M of random measure- 
ments required, in a sensing basis incoherent with the sparsity ba- 
sis, roughly scales with the sparsity, well below the number of sam- 
ples iV 2 of the signa l at Nyquist-Shannon rate: M oc K <C TV 2 
dCandes et al1l2006bl : ICandesI |2006|. 120081) . In the case of radio- 
interferometric data, the measurements are assumed to be random 
Fourier samples which represent a good sensing procedure for a 
string signal x exhibiting sparse magnitude of the gradient Vi. 
In this case, in the presence of i.i.d. Gaussian noise, the so-called 
Total V ariation minimization problem (TV e ) applies dCandes et al .1 
l2006ah . which corresponds to the minimization of the TV norm of 



4.2 Signal prior 

In a Bayesian perspective, the minimization of the TV norm of the 
signal amounts to assuming that the signal gradient follows a Lapla- 
cian distribution with no spatial correlations, and to minimizing the 
negative logarithm of this prior distribution, i.e. to maximizing the 
distribution, subject to the \ 2 constraint. This promotes compress- 
ibility in a generic sense. A Laplacian distribution is indeed more 
peaked and more heavy-tailed than a Gaussian prior for example. 
However, with realistic simulations of the string signal at hand one 
can study more precisely the signal distribution. In that perspec- 
tive, the signal gradient can be modeled at each sample point by 
the positive tail of a Generalized Gaussian Distribution (GGD): 



V{Vx\p) = 



2psT (q- 



(9) 



where V is the Gamma function, and where q and ps are respec- 
tively called shape and scale parameters. The shape parameter q 
can be considered as a continuous measure of the compressibility 
of the underlying distribution. Setting q — 2 recovers a Gaussian 
distribution, and q — 1 recovers a Laplacian distribution. Letting q 
approach yields compressible distributions. This parameter deter- 
mines the kurtosis . the ratio of the fourth central moment 
to the square of the variance (second central moment), by 



T(5q- 



[r(3g-i)r 



(10) 
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By nature this kurtosis is independent of the string tension p. The 
scale parameter ps is linearly proportional to the standard deviation 
it' Vi '(/)) of the distribution. The corresponding variance reads as: 



(Vx) 



(p) 



,T(3q- 



(ii) 



The shape and scale parameters are estimated by a moment 
method on the aggregate of the magnitude of the gradient of the 
64 test simulations. The values with two significant figures are 
q = 0.42 and s = 2.3 x 10 5 . The shape parameter value shows 
that the signal gradient exhibits a much more compressible gradient 
than a generic Laplacian distribution. Figure [4] shows the positive 
tail of the modeled prior GGD distribution V(Vx\p) for a string 
tension p = 2 x 1CP 7 . The GGD is superimposed on the his- 
togram of the aggregate of the magnitude of the gradient of the 64 
test simulations. Qualitatively, we see that the prior distribution is 
well modeled by a GGD. 



4.3 STV e problem 

Accounting for the specific statistical prior information in the reg- 
ularization of the inverse problem should enhance the reconstruc- 
tion performance, beyond the standard TV e minimization. In the 
Bayesian perspective discussed, we define a Statistical Total Vari- 
ation STV e minimization problem that consists in minimizing the 
logarithm of the GGD prior distribution modeled, subject to the \ 2 
constraint: 



mm 



MItv subject to x (y,<*>,x) 



(12) 



Again, we set e 2 to be the 99 th percentile of the \ 2 distribution. Be- 
cause q < 1 the problem is non-convex in nature, but a re-weighted 
scheme allows convergence of the algorithm thr ough iterations of 
simple TV e minimizations l lCandes et alJl2006H) . 

Note that no prior knowledge of the string tension p is for- 
mally required for the reconstruction. However we must acknowl- 
edge the fact that convergence appears to be faster when some prior 
estimation of p is used to renormalize x in the minimization term of 
relations {8} and \12\ , as expected from relation {5). As discussed 
in Section[5] such an estimation of p will also define a threshold be- 
low which one should not expect good reconstruction of the string 
network. Our estimation procedure is described below. 



4.4 String tension estimation 

By Bayes' theorem, the posterior probability distribution for p 
given the measured visibilities Q(p\y) is simply obtained from the 
likelihood C(y\p) and the prior probability distribution V(p) on p. 
For complete consistency, the likelihood should be calculated using 
the model established for the signal gradient and noise. However, 
while this model by construction accounts for the non-Gaussianity, 
i.e. sparsit y of the string signal, it ignores the spatial correlations. 

As in lHammond et all ( 120091) . we have observed that a likeli- 
hood yielding a precise localization of the string tension value can 
actually be obtained using a power spectral model. Such a model 
assumes that both the string signal and the noise arise from statisti- 
cally isotropic Gaussian random processes, such that their Fourier 
coefficients are independent Gaussian variables. Under this model, 
as the string signal and noise are independent, the observed signal 
A(x + g) has a power spectrum, as modified by the primary beam, 
of the form: 



p(A(*+ 9 )) ^ ^ = p (Ag) + p (Ax) ^ ^ _ (B) 

with P (A9) (k) defined in Section [331 and similarly P (Ax) (k, p) = 
\A(k)\ 2 p( x \k, p) . In this setting the likelihood can be computed 
easily in terms of the measured visibilities y as: 



M/2 



(14) 



The posterior probability distribution for p given the measured vis- 
ibilities thus reads as 

Q(p\y) = D- 1 V(p)£(y\p), (15) 

with normalization D = J V(p)C(y\p)dp. We take the prior V(p) 
to be flat in an interval p 6 [0, p (max ^], with an upper bound p( ml> *> 
large relative to the upper bound l[4j. 

We define an estimation p of the string tension as the expecta- 
tion value of this posterior probability distribution: 

p = E[Q(p\y)]. (16) 



5 ANALYSES AND RESULTS 

In this section, we firstly describe our numerical analyses, including 
the performance criteria of reconstruction, and discuss their results. 
We also compare the corresponding eye visibility thresholds on the 
string tension with the detectability threshold obtained on the basis 
of the PSM. We finally propose a small discussion of our results. 



5.1 Analyses 

The overall reconstruction is effective for mapping the string net- 
work if the magnitude of gradient of the reconstructed signal 
closely resembles the magnitude of gradient of the true string sig- 
nal. A simple qualitative measure of the performance is given by 
whether the string network is visible in the magnitude of gradient 
of the reconstructed signal. We define the eye visibility threshold as 
the minimum string tension around which the overall denoising and 
mapping by the magnitude of the gradient begins to exhibit string 
features visible by eye. We will augment this qualitative assessment 
of the reconstruction performance with three quantitative measures, 
namely the signal-to-noise ratio, the correlation coefficient and the 
kurtosis of the magnitude of the gradient of the reconstructed string 
signal re-multiplied by the primary beam. The kurtosis is known 
to be a good statistic f or discriminating betw een models with and 
without cosmic strings dMoessner et al.lll994l) . 

The signal-to-noise ratio is defined in terms of the original and 
reconstructed string signal gradients as: 



SNR 



(VAi,VAi) 



= -201o gl 



r (VAx — VAx) 



(17) 



where a {^Ax-vAx) and a (VAx) reS p ect j ve i v stan( j f or me stan- 
dard deviations of the discrepancy signal X7Ax — V Ax and of the 
original signal SI Ax. The standard deviations are estimated from 
the sample variances on the basis of the signal realizations con- 
cerned. With this definition, the SNR is measured in decibels (dB). 
We will consider that the reconstruction is meaningful in terms of 
signal-to-noise ratio for the values of p where this statistic is posi- 
tive. 

The correlation coefficient is also defined in terms of the orig- 
inal and reconstructed string signal gradients as: 
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Noise 


Coverage (per cent) 


FSM detectabmty 


Eye visibility 


PA - IN 


25 


7.8 X 10~ 10 


1.6 X 10" 9 


PA - IN 


50 


7.1 x lO" 10 


1.0 X 10" 9 


SA - tSZ 


25 


2.9 x 10~ 7 


6.3 x KT 7 


SA - tSZ 


50 


2.4 x 10" 7 


4.0 x 10~ 7 


SA + tSZ 


25 


5.0 x 10" 7 


1.0 x 10~ 6 


SA + tSZ 


50 


4.4 x 10~ 7 


6.3 x KT 7 



Table 1. PSM detectability thresholds and eye visibility thresholds on the 
string tension, for each of the noise and coverage conditions considered. All 
values are given with two significant figures. 



„(VAi,VAi) 



(V/li,VAi) 



(18) 



a (VAx) a (VAx)> 

where cov^ Ax ' VAx ^ stands for the covariance between X7Ax and 
S7Ax, and <r' VA:c ' and cr( VAa: ' respectively stand for the standard 
deviations of the signals VAx and VAx. The covariance and vari- 
ances are again estimated from the sample covariance and variances 
on the basis of the signal realizations concerned. We will consider 
that the reconstruction is meaningful in terms of correlation coef- 
ficient for the values of p where this statistic is also positive. Note 
that the variances and covariances discussed obviously depend on 
the string tension, but not their dimensionless ratios. For this reason 
we have omitted this dependence in the definitions l |17| > and 

Analogously, kurtoses of the magnitude of the gradient of the 
reconstructed string signal re-multiplied by the primary beam are 
estimated from the sample kurtoses on the basis of the signal real- 
izations concerned. The kurtosis of the magnitude of gradient of a 
string signal as expected on the basis of our test string signal simu- 
lation takes the value k' VA:e ' = 51 with two significant numbers. 

We compare the reconstruction results of the simple INVERT 
procedure and of the standard CLEAN algorithm with those of 
the TV E and STV E algorithms. Reconstructions are performed 
for string tensions equi-spaced in logarithmic scaling in the range 
log 10 p 6 [—9, —5], corresponding to ratio values for p of 1.0, 1.6, 
2.5, 4.0, and 6.3 in each order of magnitude. For each string ten- 
sion, each algorithm, each noise condition (PA — IN, SA — tSZ, 
or SA + tSZ) and each coverage condition (25 or 50 per cent) anal- 
ysed, we consider that the quantitative measures described above 
indicate effective performance when this performance is significant 
over the 30 corresponding simulations. 



5.2 Results 

In all noise and coverage conditions, the simple INVERT proce- 
dure is already effective at enhancing the visibility of strings by eye 
thanks to the multiplication by W 2 in Fourier space, but the result- 
ing maps are dirty. The string tensions representing the eye visibil- 
ity thresholds are listed in Table Q] These thresholds are well be- 
low the experimental upper bound in the noise condition PA — IN, 
while for the noise conditions S A — tSZ and SA + tSZ they are re- 
spectively of the same order and larger. This evolution is naturally 
due to the larger amount of noise at high spatial frequencies, where 
the signal gradient lives, in the presence of secondary anisotropics. 

The eye visibility thresholds for CLEAN, TV E , and STV E 
are the same as for INVERT. However it clearly appears that 
these more evolved algorithms offer much better reconstructions. 
This can already be acknowledged by simple eye inspection, as il- 
lustrated in Figure [5] for the magnitude of the gradient of the re- 
constructed signals re-multiplied by the primary beam in the noise 
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Figure 5. Magnitude of the gradient of reconstructed string signals re- 
multiplied by the primary beam, for a tension p = 1 X 10~ 7 , in the noise 
condition PA — IN and for a coverage at 25 per cent. The reconstructions 
considered are those associated with the simple INVERT procedure (top 
left panel), CLEAN (top right panel), TV E (bottom left panel), and STV E 
(bottom right panel). 



condition PA — IN at 25 per cent and for a string tension p — 
1 x 10~ 7 . This qualitative assessment is confirmed by our quanti- 
tative measures for the reconstruction performance in all noise and 
coverage conditions and for all string tensions. Figure |6]represents 
the SNR, correlation coefficient, and kurtosis of the magnitude of 
the gradient of reconstructed string signals re-multiplied by the pri- 
mary beam as functions of the string tension in all noise conditions 
and for a coverage at 25 per cent, for CLEAN, TV E , and STV E . 
The corresponding curves for INVERT are actually not shown for 
readability reasons as they are very far from optimal values. In this 
context, the major interest of the analysis actually resides in the 
comparison of the standard CLEAN algorithm with the TV E and 
STV E algorithms. 

Firstly, the reconstructions illustrated in Figure [5] seem very 
similar for CLEAN, TV E , and STV £ . Note however that the mag- 
nitude of the gradient of the original signal (not shown on the fig- 
ure) re-multiplied by the primary beam lies roughly in the range 
[0,3.5]. The range of amplitudes reported on the maps therefore 
indicate that the dynamic range of the reconstructions is better for 
TV E than for CLEAN, and even better for STV E . This is a generic 
behaviour that can be observed in all noise and coverage conditions 
and at all string tensions. 

Secondly, as shown in Figure [6] in all noise conditions and 
for a coverage at 25 per cent, the eye visibility thresholds identify 
with the string tension where the SNR and correlation coefficient 
get down to very small values, and the kurtosis exhibits a large 
variability over the 30 simulations around the expected value. This 
allows us to identify rigorously the eye visibility threshold as the 
string tension below which the algorithms are not effective anymore 
at recovering the string network. 

Thirdly, as also shown in Figure|6]in all noise conditions and 
for a coverage at 25 per cent, the SNR and correlation coefficient 
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Figure 6. SNR in decibels (dB) (left panels), correlation coefficient (center panels), and kurtosis (right panels) of the magnitude of the gradient of reconstructed 
string signals re-multiplied by the primary beam as functions of the string tension in logarithmic scaling in the range log 10 p S [—9, —5], for a coverage at 25 
per cent in the noise conditions PA — IN (top panels), SA — tSZ (middle panels), and SA + tSZ (bottom panels). The reconstructions considered are those 
associated with CLEAN, TV e , and STV e . The vertical lines on the curves represent the variability at one standard deviation of the estimated statistic across 
the 30 test simulations considered (these lines are not visible where smaller than the width of the curves). The black dot-dashed vertical lines represent the eye 
visibility thresholds. The black dot-dashed horizontal lines identify either the limit of zero SNR, zero correlation coefficient, or the kurtosis of the magnitude 
of gradient of a string signal as expected from the test string simulation itself: k^ Ax ^ = 51. 



are always significantly larger for TV e than for CLEAN. They are 
always at least as good for STV e as for TV E and often signifi- 
cantly larger. However, it appears that the kurtosis is less optimally 
recovered with TV e than with CLEAN. But the STV e approach 
regularizes the kurtosis significantly better than CLEAN in most 
cases. 

Note that these last two conclusions also hold for a 50 per cent 
coverage, hence suggesting their stability and robustness relative 
to the coverage. The reconstruction performance is always signifi- 
cantly better at 50 per cent than at 25 per cent, but the improvement 
is not drastic. This suggests that a small number of visibilities can 
already lead to good reconstructions of the string signal. 

5.3 PSM detectability threshold 

In the perspective of the evaluation of the performance of the recon- 
struction algorithms, it is natural to ask whether the requirement for 
a complete reconstruction of the signal strongly limits the lowest 
string tension accessible in comparison with the lowest detectable 
string tension on the basis of the PSM. 



In this regard, for any possible string tension p, the probabil- 
ity distribution for the estimate p on the basis of the PSM may be 
obtained from simulations. One can identify the critical value po 
such that, for null string tension, the probability that the estimate p 
is larger than po is small: p(p ^ po) — a, for some suitable posi- 
tive value a much smaller than unity. The test for the hypothesis of 
null string tension is then defined as follows. For estimated values 
p ^ po, the hypothesis of null string tension may be rejected with a 
significance level a. On the contrary for estimated values p < po, 
the hypothesis of null string tension may not be rejected. 

We define the PSM detectability threshold p* such that, for 
the string tension p* , the probability that the estimate p is larger 
than po is large: p(p po) — 1 — j3, for some suitable positive 
value p much smaller than unity. Consequently, for string tensions 
larger than p* , the probability of rejecting a null string tension on 
the basis of the hypothesis test defined is larger than 1 — p, called 
the power of the test. The value p* is the smallest string tension that 
can be discriminated from the hypothesis of null string tension for 
given values of significance and power of the test. It may thus be 
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understood as a detectability threshold determined on the basis of 
the PSM. 

The PSM detectability thresholds in the various noise and cov- 
erage conditions considered are reported in Table Q] for a ~ (3 ~ 
0.01. In the noise condition PA — IN, detectability thresholds are 
well below the experimental upper bound, while for the noise con- 
ditions SA — tSZ and SA + tSZ they are respectively of the same 
order and larger. In all cases, the PSM detectability threshold p* 
is only slightly below the eye visibility threshold associated with 
the reconstruction algorithms discussed, by a factor smaller than 2. 
This discrepancy between these two thresholds confirms that the 
detection problem alone can be solved with the PSM only down to 
slightly lower string tensions than th e more difficult reconstr uction 
problem. As already emphasized in lHammond etafl ( l2009h . it is 
one thing to estimate a single global parameter such as the string 
tension on the basis of a PSM, but quite another to explicitly recon- 
struct the string network itself. These comparison results are thus 
already very positive. 

5.4 Discussion 

Our results strongly support the idea of the superiority of the com- 
pressed sensing approaches for a string signal reconstruction from 
radio-interferometric data, in particular when specific statistical 
prior information is introduced in the global minimization problem 
defined for reconstruction. However, at the high angular resolution 
considered and in the noise condition SA + tSZ, the eye visibility 
thresholds, as well as the PSM detectability thresholds, are above 
the current best experimental upper bound on the string tension. 
One first approach to circumvent this difficulty in the perspective 
of the reconstruction of a possible string network from forthcom- 
ing radio-interferometric observations is to consider an observation 
frequency around 217 GHz so as to eliminate the thermal SZ effect 
and to fall in the noise condition SA — tSZ. But even in this case 
the eye visibility thresholds remain slightly higher than the experi- 
mental upper bound. 

In this context, recent work by the authors might help to dras- 
tically enhance the signal reconstructi on performance in t he con- 
text of compressed sensing approaches dWiaux et alj2009bh . Radio 
interferometers with small field of view and baselines with non- 
negligible and constant component in the pointing direction should 
be considered. In this context, the visibilities measured essentially 
identify with a noisy and incomplete Fourier coverage of the prod- 
uct of the planar signal multiplied by the primary beam with a linear 
chirp modulation. Theoretical and numerical results show the uni- 
versality of the corresponding spread spectrum phenomenon rela- 
tive to the sparsity basis, in terms of the achievable quality of re- 
construction through global minimization problems such as TV e 
and STV e . Preliminary numerical results regarding the application 
of such a framework to the reconstruction of a string signal in the 
CMB look very promising. 

In another line of thoughts, global component separation tech- 
niques might be envisaged in order to simultaneously extract all 
non-Gaussian components of the CMB temperature data, includ- 
ing the SZ effects and the string signal. 

Finally, let us recall that the string signal simulations used for 
modelling the prior statistical distribution of the signal gradient as 
GGD's and for computing its power spectrum are realistic. Never- 
theless, further analysis should be performed in order to study the 
stability of the STV S reconstruction quality relative to a discrep- 
ancy between the statistical distribution of the true signal and our 
model. Analyzing the difference between the TV e and STV e re- 



construction qualities already provides useful information in that 
respect, as the TV e minimization problem captures the essence 
of the compressiblity of the signal gradient, but with a non-exact 
GGD shape parameter q = 1. In this context, one can argue in 
favour of a relative stability of the reconstruction quality in terms 
of the SNR and correlation coefficient as these measures are indeed 
smaller for TV e than for STV E , but still larger than for CLEAN. 
This stability is lost in terms of the kurtosis though. 



6 CONCLUSION 

The STV e algorithm proposed for the reconstruction of a string 
signal from arcminute resolution radio-interferometric data of the 
CMB is designed in the framework of compressed sensing. It no- 
tably relies on a model of the prior statistical distribution of the 
signal, fitted on the basis of realistic simulations. The algorithm 
shows superior performance relative to the standard CLEAN al- 
gorithm. Secondary CMB anisotropies strongly hamper the signal 
reconstruction though, and further work is still required in order to 
obtain good performance at string tensions below the current ex- 
perimental upper bound. 
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